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In this paper we examine the validity of the linear perturbation theory near a bounce in the 
covariant analysis. Some linearity parameters are defined to set up conditions for a linear theory. 
Linear evolution of density perturbation and gravitational waves have been computed previously. 
We have calculated the vector and scalar induced parts of the shear tensor. For radiationlike and 
dustlike single fluid dominated collapsing Friedmann-Lemaitre-Robertson- Walker background it is 
shown that the linearity conditions are not satisfied near a bounce. 

PACS numbers: 04.20.-q, 98.80.-k, 98.80.Jk, 98.80. Bp 



I. INTRODUCTION 

The cosmological models for a collapsing universe 
approaching an expanding one through a nonsingular 
bounce or its cyclic repetitions are being studied over 
a long time For example, in a dustlike scalar field 

dominated contracting Friedmann-Lemaitre-Robertson- 
Walker (FLRW) background both the curvature pertur- 
bation £ and the Bardeen potential $ have been observed 
to grow, whereas in an expanding phase £ remains con- 
stant and $ decays Q. Such growing modes of the per- 
turbations raise doubts on the validity of the linear per- 
turbation theory close to the bounce [l(| ■ However, some 
recent papers [ll[ , [l2| claim that there exists a choice of 
gauge in which the linear approximations remains valid, 
whereas in other commonly defined gauges such approx- 
imations become invalid. 

The power spectra of 4> and C are given in 12] for a 
contracting phase (just before the bounce) and for an ex- 
panding phase (long after the bounce) . $ has a constant 
mode and a growing mode close to the bounce. For ex- 
ample, in a dust dominated contracting flat FLRW back- 
ground, $ behaves as 



$k(r)) = Ci(k) + C 2 (k) V - 



(1) 



where r\ is the conformal time and a is the scale factor, 
a oc {—if) 2 . Due to the growing mode, $ becomes very 
large near the bounce. On the other hand, although £ 
grows near the bounce, its amplitude remains always less 
than the constant mode of <£>. If the perturbations are as- 
sumed to be small at some initial time far away from the 
bounce then the constant mode of $ and hence £ remain 
small near the bounce. In order to examine the effect of 
these growing modes of $ and £ on the validity of linear 
approximations, one needs to compare metric perturba- 
tions with their background values. Metric perturbation 
about an FLRW background is written as 



(2) 



Let us consider the scalar perturbation cf>, defined by 
goo = — 1 + 20. A necessary condition for the pertur- 
bation to be linear is that <h <C 1. In the Newtonian 



gauge, (j> is equal to the Bardeen potential Now from 
the behavior of $ it is evident that <f> ^> 1, close to the 
bounce. On the other hand, in the uniform curvature 
gauge 0=!(l + ui)£ina single fluid dominated phase. 
Since £ is small close to the bounce, 4> also remains small. 
It shows that the perturbation remains linear close to the 
bounce if one uses the uniform curvature gauge instead 
of the Newtonian gauge. Similarly, at the bounce, where 
Hubble's constant goes to zero, neither the Newtonian 
nor the uniform curvature gauge but the synchronous 
gauge preserves the linearity of the perturbation. This 
shows that there exists at least one gauge in which the 
perturbations remain small near bounce so that linear 
perturbation theory holds. But in other gauges pertur- 
bations grow. 

The purpose of this paper is to analyze the problem 
in a fully covariant approach first proposed by Hawking 
(l3j and further developed in [lij]. In this approach, the 
dynamical variables are fully gauge invariant and hence, 
the analysis is completely independent of the choice of 
gauge conditions. Moreover, in contrast with the per- 
turbative expansion of Einstein equations used in stan- 
dard perturbation theory, the dynamical equations are 
exact and nonlinear. One can relate the variables used 
in the covariant analysis with the gauge invariant vari- 
ables used in the standard perturbation theory. For ex- 
ample, the vector £„, defined in Eq. ([26| of Sec. IIII1 can 
be seen as a generalization of the curvature perturbation 
£, used in linear theory. In particular £ a coincides with 
the usual £ for long wavelengths, but the two quantities 
differ on small scales, where spatial derivatives cannot be 
neglected. 

We consider collapsing radiation and dust dominated 
FLRW backgrounds. First, the linear perturbation equa- 
tions are solved in these backgrounds near the bounce. 
Then these solutions are used to compare the linear and 
nonlinear terms in the full nonlinear perturbation equa- 
tions to investigate the validity of the linearized approx- 
imations. Our discussion is general and not specific to a 
particular model of bounce. 

So far it has been checked whether perturbations are 
sufficiently small compared with their background quan- 
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tities. If that criterion holds, second order perturbations 
are assumed to be even smaller and linear perturbation 
theory is considered to be a good approximation. But 
this test will not work in covariant formalism as the 
background values of the gauge invariant perturbations 
are zero. In the present work, we have checked whether 
higher order perturbations are truly smaller compared 
with 1st order ones. We will show that the answer is not 
positive for all perturbation modes near bounce. 

The outline of this paper is as follows. In Sec. UTJ 
we set up our notations and give a brief introduction 
to the covariant perturbation theory. In Sec. IIIII we list 
all the gauge invariant variables. In Sec. IIVI we give the 
full nonlinear perturbation equations and their linearized 
versions close to a flat FLRW background. Also the con- 
ditions for validity of the linear approximation are set up. 
In Sec. [V] we present the solutions of the linearized equa- 
tions for a radiation and dust dominated background and 
finally, in Sec. I VII we compare the first and higher order 
terms in full nonlinear equations. 



II. THEORY OF COVARIANT 
COSMOLOGICAL PERTURBATION 

In the standard perturbation theory [l5j],[l||, [l?} we 
consider an idealized model of universe S, usually a 
FLRW one with metric g ab , energy density //, pressure 
density p, etc. Then we perturb the model to obtain a 
realistic universe S with corresponding quantities g ab , \± 
and p. The perturbations at each space-time point q is 
defined as the difference of these quantities at q. 

Sgab = 9ab ~ 9~ab, Sfi = fl~ /2, Sp = p - p. (3) 

But S and S have different manifold structure. Mathe- 
matically these means we have to define a map $ from S 
to S. Then we can pull back each quantity 11 in S to its 
image in S. 

= /*(#(«)), qeS, q = 4>(q)eS. (4) 

Then the perturbation is defined as 

5»(q) = - fi(q). (5) 

There are coordinate freedoms in both the manifolds S 
and S and these choices determine the map $. These 
arbitrariness in the $ makes the definition of $*/i(q) and 
hence Sfi ambiguous. This is the gauge freedom in /i. 
To study the evolution of perturbations we have to fix 
the gauge i.e. specify the map then there will be no 
ambiguity in 5[i. 

But, S is the realistic universe and S is our theoretical 
idealization. Given a physical universe S we can not 
construct the background model S uniquely without any 
further information. 

In [l3j], [l4| an alternative representation to stan- 
dard perturbation theory is given. This representation 



involves quantities, whose background values are zero. 
Suppose / is such a quantity, defined on S. Then we 
can consider / itself as the perturbation variable and it 
is completely gauge invariant. Consequently, the evolu- 
tion equations for / are fully covariant. The quantities, 
which do not vanish on background space-time, are called 
zeroth order variables. By the statement, "The back- 
ground space-time is flat FLRW," we mean that zeroth 
order quantities are governed by the Friedmann equations 
with zero 3-curvature. 

Now we consider a general perfect fluid flow in curved 
space-time [l8j,!lS]. In cosmological context, there will 
always be a preferred family of world lines representing 
the motion of typical observers. Let the four velocity of 
those fundamental observers be u a . 

dr a 

u a = ^-, u a u a = -l, (6) 

(XT 

where r is the proper time along the fundamental world 
lines. The metric of the 3-space, orthogonal to u a , is 

Kb = 9ab + u a u b =$> h\h h c = h a c , h ab u b = 0. (7) 

Time derivative of any tensor S a b " along fluid flow lines 
is defined as 

.V'; ;. V, .S"-;;; . (8) 

The first covariant derivative of four velocity can be writ- 
ten as 

Vfetto = -0h a b + (T a b + U> a b - U b U a , (9) 

Where the trace V a u a = is the expansion, traceless 
symmetric part, a a b = C( a 6) is t ne shear tensor, anti- 
symmetric part, uj ab = uj[ab\ is the vorticity and u a is the 
acceleration. We can define a scale a(r) along each world 
line as 

= 3-, (10) 
a 

a can be determined up to a multiplicative constant. 

<7 ab U h = = UJ ab U b (11) 

The space-time curvature (Riemann) tensor R a bcd, de- 
fined by 2Vr a VMt! c = R a bcdV d is made up of Ricci tensor 
R a b and the tracefree Weyl tensor components C a bcd- 

Rabcd = Cabcd + 9a[c.Rd]b ~ 9b[c.Rd]a ~ -^9a[c9d]bR (12) 

The trace part R a b is determined by the Einstein equa- 
tion 

Rab ~ ^Rgab = nT ab - (13) 

k = 8ttG is the Plank length squared. The energy mo- 
mentum tensor T a b for a perfect fluid can be written in 
terms of energy density /i, pressure density p and u a , 

Tab = (M + P)u a Ub + P9ab = f-lUaUb + ph ab . (14) 

The energy-momentum tensor must satisfy the conserva- 
tion equation, 

V b T ab = 0. (15) 
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III. GAUGE INVARIANT VARIABLES 

A FLRW model is a per fect fluid space-time character- 
ized by the conditions [14j . 

<Jab = U ab = U a = 0, (16) 

which implies 

/i = //(*), P = P(t), = 0(t), (17) 
where t is the cosmic time, defined by u a — — V tt 

t; and, 

Weyl tensor vanishes 

C abcd = 0. (18) 

So this space-time is conformally flat. 

From the above characterization, we can list some sim- 
ple gauge invariant quantities (20j : 

(1) Shear, vorticity and acceleration, 

°ab = (h {a h d - h ab h cd )W d u c , (19) 

uj ab = h£h b fV d Uc, (20) 

u a = u b V b u a . (21) 

Rotation vector is defined as uj a = \^abc^ bc '■ tabc is Levi- 
Civita tensor in 3-hypersurface defined by e abc — r] abcd u d . 

(2) "Electric" and "magnetic" parts of the Weyl tensor, 

Eab = C acbd u c u d , Hab = \c acpq rf\ d u c u d . (22) 

(3) Spatial gradients of the energy density, pressure den- 
sity and expansion, 

X a = Kh a b V b fi, Y a = nh b V b p, Z a = h b V b 9.{23) 

There are many other gauge invariant quantities, which 
vanish in the background model. We should mention two 
of them those will appear in the evolution equations. 

A = V a ii Q , A a = h b V b A. (24) 

It is evident that fi, p and 9 are zeroth order variables, 
whereas those defined in (rj~9"l) - (|2"4")) are called first order 
variables. We can construct higher order variables from 
first order ones as 

. ,2 1 , , , ,ab 2 1 „ ab „ TP&b „ Tjab (r.r\ 

w = ^abOJ , er = -cr ab <7 , a ab h , cr ab H (25) 

The comoving fractional density gradient, defined as 
T> a = a-^-, is similar to the variable 5 = 5^//i, used in 
standard perturbation theory and is observable in prin- 
ciple. So we will concentrate on the equations of X a . 
The gauge invariant quantities defined in (|19 j) -(|23 [) form 
a closed set of equations, as will be seen in the next sec- 
tion. 

To demonstrate the relation between covariant and co- 
ordinate based approaches let us define another variable 

Ca as [m, m 

(a = W a + - Xa , W a = h b V b log a. (26) 
3k(^ + p) 

A cosmological space-time, close to the flat FLRW geom- 
etry can be described by the metric in comoving coordi- 
nates, 

ds 2 = -(1 + 2<t>)dt 2 + 2a(a tJ B + S i )dtdx i 
+ a 2 ((l + 2i/;)% + 2d l d J E + 2d (l F 3) + h ij )dx i dx j {27) 



The matter content is a perfect fluid with energy den- 
sity /j,(t 7 x l ) = p,(t) + Sfi(t, x l ), pressure p(t, x l ) — p{t) + 
Sp(t, x l ) and 4- velocity u M = + 8u^ . 

Suq = —(/>, Siii — di5u + Suj- (28) 

At linear order, spatial components of £ a are 

C^C S , C S = --^. (29) 
a fi 

Up to linear order 8a/ a can be written as 

= $ + -W 2 E + - f dtW 2 Su. (30) 
a 3 3J 

So, C S coincides with familiar definition of comoving cur- 
vature perturbation £ = ijj — when the spatial gra- 

dients are negligible. For adiabatic perturbations, £ s is 
conserved at all scales whereas £ is conserved only for 
large wavelength modes. 



IV. DYNAMIC EQUATIONS AND 
CONSTRAINTS 

A. Exact equations 

The time and space components of the energy- 
momentum conservation equation (|15[) leads to 

/i + 0(/i + p) = 0, (31) 
K(fi + p)u a + Y a = 0. (32) 

The Raychaudhury equation gives the evolution of 9 
along fluid flow lines Q , Q 

9 + )-9 2 - A + -k((i + 3p) + 2(cr 2 - uj 2 ) = 0. (33) 

To solve the zero order background equations (j3Tj) and 
(|33[) we need the matter equation of state in the form 
p = p(n). Then (|23|) and (|32|) imply that u a and Y a can 
be written in terms of X a . 

The equations for X a and Z a are given in [ljj 

a- 4 h b (a 4 X b ) = -«0i + P)Z a - (oj\ + a\)X b , (34) 
a- 3 h b (a 3 Z b ) = Kii a - ^X a + A a + 2h b V b {u 2 - a 2 ) 

-{uj b a + a b a )Z b , (35) 

where, 

K = Kfi- \9 2 + A + 2(uj 2 - a 2 ). (36) 

Equations for four other variables ui ab , a ab , E ab and H ab 
are given in [l3| 
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a- 2 h a c h b d (a 2 Lj cd ) 
a- 2 h a c h b d (a 2 a cd ) 
a- 3 h a c h b d (a 3 E cd ) 



hfh h d \/ [d ii c] 



-E ab 



2cr c [ a w 



(bU a ) 



b] 



LU ac LU 



+ ^h a b(cr 2 - w 2 ) + U a U bl 



-CurLff ob - -k{h + p)(T a b + E c {a uj b)c + E c (a a b)c + e acd e bpq a cp E dq - 2H c {a e b)cd ii d , 



i h a c h b d (a 3 H ed ) = CurlEah + H c {a uj b)c + H c {a a b)c + e acd e bpq 



a cp H dq _ 2H . 



\a e b)cdU 



(37) 
(38) 

(39) 
(40) 



We have used following notations [25| : 

\ab) — 2^ ab + ^ ba )' \°-b] = l^i^ab — \a) , 

•"V&) = h^h b d {\^ cd ) - -h cd X e e ), 
Cm\X ab = h e {a e b)cd \7 d \ e c . 

There are also constraint relations which must be sat- 
isfied at some initial time on each world line 1131. 



h a c V b {uj b c + a\) ~ u b {u ab + a ab ) = 


\z a , 


(41) 


V a w a = 2u a oj a , 




(42) 


Curlw ab + Curkr ab = -H abl 




(43) 


h a c V b E bc + 3H ab u; b — e ab cP b d H cd = 


(44) 


h a c V b H bc — 3E ab u b — e abc a b d H cd = 


K(p+p)0J a 


(45) 



B. Linearization about FLRW model 

Let us consider the universe to be almost FLRW and 
linearize the perturbation equation. Since /i, p, 9 are ze- 
roth order variables, we consider only zeroth order terms 
of (|3Tj) and (f33|) . Then the Raychaudhury equation be- 
comes 



1 



1 



+ -e 2 + -k(jx + 3 P ) = o. 



(46) 



Now a can be interpreted as the Robertson- Walker scale 
factor. Linearized forms of the perturbation equations 
become 



a- 4 (a 4 X a 
a- 3 (a 3 Z a 

-2/ 2 

a (a Lu ab 

-2/ 2 

a [a (Jab 
a~ 3 {a 3 E ab 
3 (a 3 H ab 



K(p+p)Z a , 

\-A a , 



1 

~2 Xa 



= -E a 



1 



= -CurLff afc - -n(p +p)a ab , 
= CurLB a6 . 



(47) 

(48) 

(49) 
(50) 

(51) 

(52) 



Since background FLRW universe is flat, 1Z is a first order 
quantity. That is why the 1st term in the right hand side 
of (l3"5t is neglected in ([4*5]). Also we have the constraints, 



W b {u b a + a b a ) 



V a w a = 0, 



(53) 



Uab = -Curlw ab - Curl<7 a6 , 
1 



V h E ab 



V H ab = n(ii+p)uj a . 



(54) 
(55) 



C. Condition for linearity 

In the standard perturbation theory, linearization is 
justified if the perturbations are small with respect to 
corresponding background quantities. The gauge invari- 
ant variables defined in Sec. Mil are nonlinear and so far 
we have not assumed any "smallness." However when we 
consider the linear evolution, some conditions have to be 
imposed. They cannot be compared with corresponding 
background quantities, because those are zero from the 
definition of gauge invariant variables. The definition of 
dimensionless variable [26j is not unique, because one can 
multiply any power of scale factor a, and that variable 
remains dimensionless. The natural way is to demand 
that higher order terms of perturbation equations remain 
small with respect to first order terms. 

Let us concentrate on the evolution of density pertur- 
bations X a . It is observed that X a is coupled with Z a 
up to linear order in Eqs. (|3"4")l and ([331) . To compare 
nonlinear terms of those equations with linear terms we 
define following parameters: 



£i 



£3 



£e 



£3 



VaXb\ 



\K(fJ, + p)Z a \' 



£2 = 



W\Xb\ 



K([J,+p)Z a \ ' 



£6 = 



\KUg\ 

\l X "\ ' 

\<z b \ 
I \ X A 

\KUg\ 

\Aa\ ' 
\" b gZb\ 

\AJ 



£4 



£7 



£4 



\2h b V b w 2 \ 

\l X a\ 
V\Zb\ 

I 1 X I ' 
l2 A a| 

\2h b V b W 2 \ 
I A, I 



£5 



\2h b V b a 2 

\i x <A 



\2h b V b a 2 \ 



£5 



14.1 



£7 = 



W\z b \ 



\A n 



(56) 



The linear perturbation theory is valid for the solutions 
of (l34|) and (|35j) . if the following conditions are satisfied 
throughout the regime under consideration: 

(1) £i,£ 2 < 1, 

(2) e 3 , £4, £5, £6, £7 < 1, and/or £3, £4, £5, £6, £7 < 1- 
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V. SOLUTION OF LINEARIZED EQUATIONS 
FOR SPECIFIC MATTER DESCRIPTION 



Now we will concentrate on adiabatic perturbations in 
a collapsing FLRW background, for "Dust" and "Radi- 
ation" dominated cases. The solutions of coupled equa- 
tions (j47|) and (|48|) for X a and Z a are obtained inde- 
pendently. Then u a is found from (|32[) . The w ab also 
obtained from (|4"9")) . Then we can solve (f5Ul)- (|52p using 
uj a b and u a as a source. 

Let us consider equation of state of the dominating 
single fluid near bounce is 



Ya = wX a . 



(57) 



p = Wfl 

Zeroth order quantities, obtained from (j3"Tj) . (|46|) and (flO|) 

are 



M 



t 



a 3(l+«) ' 



j 3(l+to) 



1 



1 + tot 
(58) 

We assume the bounce occurs at t = —t b such that 
tf, is of the order of Plank scale. The initial time t = 
—t,,, is well inside the regime dominated by the single 
fluid considered here such that, a b = a(— t b ) « 1. The 
constants M and t* have the following relation: 



KMti 



3(1 + w) 2 ' 

The acceleration is expressed in terms of X a as 

w X a 
u = . 

1 + W KfJb 



(59) 



(60) 



Now substituting ([SU|) in (|4"5)l and using the background 
quantities from (|58|). we extract the second order differ- 
ential equation of X a from (|4"T)) and (I48[) . 



-X„ H 0X n 



ll + 3«; 
6 



2 X a - wh a b V b V c X c = 0. 

(61) 



In principle, we can solve (|6Tj) and obtain Z a from (|47l) . 
Solutions for density perturbation T> a = for radia- 
tion and dust dominated FLRW background are given in 
[27j and [H| respectively. 

The w a b is easily obtained from (|39|) . Since X a and 
hence ii a is spatial gradient of a scalar, the right-hand 
side of (|3ni) can be simplified for a torsion-free manifold. 



l + W KfJ, 

w 1 
1 + W KjJi 



lh a c h b d V [d h^ e ^ 



~ h a h b V[ d U G 



l + W (1 
W0uj ab 



Substitution of this in (|49[) results in simple evolution of 

Uab- 

oJab = n ab a- 2+3w , (62) 
where f2 ab = 0. Let us consider another important quan- 
tity, 



r a = V UJ ab 



(63) 



For any first order (0,2) tensor X a b in FLRW background, 
the time derivative and spatial divergence do not com- 
mute; rather, they satisfy the following relation: 



(h a b V c X bc ) = hy'aia^Xbc). 
Using (f62|) and (|64|) we get the evolution of r ai 
r a = R a a- 3+3w , R a = 0. 



(64) 



(65) 



The other three variables a ab , E a b and H a b are trace- 
free symmetric tensors. Equations ([50| - (|52p can be recast 
in the form, 



l-3io 



2 &ab = V( (a 2 (a 2 ?i h ))) - V( a Z t 



2V( Q r b) , 
1 



A£ ab + Is o6 + H(i - w )P Eah il + W f 8 +3w) e 2 a ab = --(1 - w)V {a X b) , 

AH ab + 7 -H ab + hi - w)e 2 H ab = -\{1 + w)0 2 V/ o um, 
3 3 6 N ' 



(66) 
(67) 
(68) 



where, AA„ 



V 2 A ab 



h a c h b d hPiV ' P {h c r h d s hjV ' t \ rs ). In deriving (|56 ]) -(P ]) . we have used following 



identities for B ab = h£h b \B cd and v a — h b v b 



(CurLB Qb ) = Curla(a 1 B ab ), (V ( a v b ) j = V( Q a(a 1 v b) ), 



(69) 



Curl 2 S ab = l -(h ab V 2 - h ( c a h b d V c V d )B e e + ^V {a h b)c V d B d 



2h {a h b d V d V e B [t 



V 2 B, 



(ab)- 



(70) 
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The homogeneous parts of the differential equations Thus we obtain an equation for Fourier mode X(k,t), 
obtained by setting X a ,Z a ,ui a equal to zero 



X(k, t) + Hx(A, t) + ( + £2 ) x ^ f ) = °- ( 74 ) 



9_ _fc 2 
2*2 + 3a 1 



give the pure tensor perturbations or gravitational waves 
[13| . |28j . However there are inhomogeneous parts of 
those equations. We have presented the solution of the 

inhomogeneous equation (j^J) for a ab , using X a , Z a and The first term in the coefficient of X in (JT^J) grows as 
uj a as source. In particular, we will consider superhori- t~ 2 — a -4 , whereas the second term grows as a~ 2 . So 
zon modes which are responsible for large scale inhomo- except for some very short wavelength modes the first 
geneities. 



A. Radiation, "w = h r 



term dominates over the second, near the bounce. 
In terms of conformal time, 

f dt 

n = — / — = — 2t. lf a, 
J a 

the general solution of ([74)) becomes 

krj 



(75) 



For radiation dominated collapsing FLRW background X(k,rj)—7] 2 ICi(fe)Y|l — = 

— I +/+ U/2 a _ q /o+l ftiTTl + +U„ f^,-™ v \ vo 



[a = (-V^*) 1/2 , = 3/2i], (EU) takes the form, 



1q + Yt ±a + 2^ X ° ~ ^> b V c A c = 0. (71) 



Since X a is constructed from a scalar /1, we can expand 
X a in Fourier modes in terms of spatial harmonics, de- 
fined in the Appendix, 



A a = ]TA(M)Qi 0) . (72) 

k 

The last term in the left-hand side of ([71]) reduces to, 



C 2 {k)j3 



kr\ 

71 



(76) 

where J is Bessel function and Y is Neumann function. 
For superhorizon modes, 



- < \H\ =*► 2fct„a< 1. 
a 

fc>; 2kt„a 



(77) 



In that limit, — ^ = << 1, we can use the asymp- 

totic expansion of Bessel functions for a small argument. 



A(M) - X^(k)a- 6 +XW(k)a- 3 . 
Then, using (j47|) . we obtain 



where, 



Z(M) = Z«(fc)a~ 4 + Z^(k)a-\ 



(78) 



(79) 



h„ b V b V c A c 



3 a 



i^A(fc,t)^ Q b V b V c /i c d V rf Q 
i£5*(M)Q< >. 



(0) 



z W(k) = le- l x^{k), zv\k) = -\e- l xW{k). 



3 * — ' a 

k 



(80) 

is the value of expansion at initial time surface t = — t*. 
The equation (1661) for radiation dominated background 
(73) takes the form, 



5 1 
'dab + -6& a b - W ab = -j-j^V {a (a- 2 (a 6 X b} )) - V \ a Z b) - 2V (a r b) . 



(81) 



o a b is a traceless symmetric tensor and one can decom- 
pose it in three parts, 

(Tab = vlb + a ab + a ab 

= E- T ^ ) +E^ ) +E^i° ) - ( 82 ) 

First two terms in the right-hand side of Eq. (|8"Tj) are 
constructed from scalars, whereas the third term is orig- 
inated from a divergenceless vector. 



V {a (a- 2 (a 6 X b) )) = Y,(a- 2 (a 6 X(k,t)))V {a Q 



(0) 
b) 



-Y^ke^ 2 \k)a- 2 Q^. (83) 



,(0) 



V (a Z b> =Y J Z{k,t)V {a Q { b) 

k 

= -J2k(z^(k)a- 5 + Z^(k)a- 2 )Q^.m 



From (|65p we can expand r a in Fourier modes for w = 3 , 
r a = £ r(A, = E ^( fc )^ 2 Qi 1} . (85) 
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V (a r b> = ^r{k,t)V {a Q$ = -^fci?(fc)a- 3 gi 1 b ) {86) 

fc k 

Using the Fourier expansion, ([82" ]) .([83 ]) . ((84)) and ([86]). the 
Eq. (|8T|) can be decomposed into three parts, 

5 fc 2 
a T (fc, t) + -6& T (k, t) + -^cr T {k, t) = 0, (87) 
3 a z 

a v (k,t) + -9& v (k,t) + ^<r v (k,t) = V(k)a- 3 , (88) 
3 a 1 

a s (k,t) + le& s (k,t) + ^a s (k,t) 
3 a z 



S {1 \k)a- 5 + S (2) (fc)a~ 2 ,(89) 



where 



V = 2fci?, S« = 5^ = -^-XW. (90) 

3a* 120* 



Using conformal time (]75[) . the propagation of gravita- 
tional waves are obtained from ([87]) [281 . 



a T (fc ;?7 ) = rj-i (X»i(*)yj (-foj) + L> 2 (fe)J f (-fcr?)) 



In the supcrhorizon limit —/cry = 2kt*a -C 1, 
a T (fc,i) = s^! ) (fc) a - 3 + S^ 2) (fc). 

.2 

In this limit we can omit the ^ term in ([87]) - 
the solutions of ([88]) and ([89]) are 

cr y (fc, f) = Ey (fc)a, 

a s (fc,t) = 4 1) (^ 1 + 4 2) W a2 ' 



(91) 

(92) 
Then 

(93) 
(94) 

Propagation of acceleration ii a and its spatial derivatives 
A, A a are obtained from ([50]). (IM1) and (175]) . 

"« = -i^ 2 E (^ (1) (fc)«- 2 + ^( 2 )(fc)a) Q<°>, (95) 
fe 

A = -^- 2 ^ fc (X( 1 )(fc)a- 3 + X( 2 )(fc)) Q(°), (96) 
^ = W E fc2 (^ (1) (fc)«- 4 + XW(k)a-A Qi°).(97) 



B. Dust, "to = 0" 

In a dust dominated background, the density pertur- 
bation is scale invariant, as evident from ([6"T]) . Also the 
acceleration is zero. The behavior of X a and Z a are [3| 

X(k,t) = X^(k) a - n / 2 + X^(k)a~ 3 , (98) 

Z(k,t) = Z^\k)a- 4 + Z^\k)a- 3 / 2 . (99) 

However, the shear is not scale invariant. For superhori- 
zon wavelengths, pure tensor components of <j ab are [28[ 

a T (k,t) = T,^\k)a- 3 + T,^\k)a- 1 ^. (100) 
The vector and scalar components of a ab are found to be 
a v (k,t) = Sy(fc)a- 1 , (101) 
CT 5 (fc,<) = El! ) (fc)a- 3 loga + E^ 2) (fc)a 1/2 . (102) 



VI. COMPARISON OF LINEAR AND 
NONLINEAR TERMS 

We have seen that each variable has a growing mode(s) 
near the bounce. But this growing mode(s) can be ab- 
sorbed by a mere redefinition of that variable. For exam- 
ple, if a variable / behaves as / ~ a~™ near the bounce, 
then we can construct a new variable / = a m f, such 
that m > ri, which remains finite. So the growing modes 
of perturbations do not rule out their linear evolution. 
In order to investigate whether the perturbations remain 
linear near the bounce one needs to compare linear and 
nonlinear terms of the perturbation equations using the 
behavior of linearity parameters defined in Sec. II V CI 



A. Radiation dominated case 

Near the bounce, we will consider only the dominating 
modes of different variables. 

X a = X a a- 6 , Z a = Z a a- 4 , (103) 
u a = U a a- 2 , A = Aa- 3 , A a = A a a- 4 , (104) 
(J ab = Q ab a~ 1 , er afc = £ ab a~ 3 . (105) 



The higher order value of curvature perturbation (]36|) is 

1Z = A + 2{ui 2 — a 2 ) 

= Aa- 3 + 2Q 2 a- 2 - 2£ 2 a~ 6 

= -Ka,- 6 , (106) 

where, X a , Z a , U a , A, A a , Q ab , S afc and H arc time 
independent. 

Using the commutation of time and spatial derivatives, 
h a b V b u 2 = fl a a- 3 , h a b V b a 2 = £ a a- 7 , (107) 

with fl a = S a — 0. Now we can calculate the linearity 
parameters defined in (J56j, 



Si 



£3 



£3 



£6 



\%KMZ a 
\KU a \ _ 
1 2-^0 1 

I \ x * 

\RUa\ 
| la | 

\n b gZb 
\AJ 



■a, e 2 



?4 



\kMZJ 



\hX a 



■ft i £5 



(108) 



L a, £7 



a 4 , £4 



-a , £ 7 



1 2fi a I 



'6 _! 
-r-ft 



-a, £ 5 



|2S a |._ s 



|^ 6 a-^fc| —J 

UJ 



(110) 

(111) 

(112) 



Let at some time slice t = —t\, such that t* 
ti ^> tj>, linearity conditions are satisfied. So £ 2 (— ii) = 

a," < 1, If we consider another time slice 



, 3 ,5MZ„ I 1 

t = — 12, which is close to t b , i. e. t\ 3> ti tb, then 



£2(-*2) = £2(-*i) 



01 

02 



(113) 
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Since 02 -C ai, the parameter £2 may become order 1 at 
t = —t 2 and the condition (1) no longer holds. Similar 
arguments can be given for £3, £5, £7, £3, £5, e~e and £7. 



B. Dust dominated case 

The evolution of variables in dust dominated collapsing 
FLRW background near the bounce, 

X a =X a a- 11 / 2 , Z a = Z a a- 4 >. (114) 
u a b = V-aba~ 2 , o" ab = £ ab a~ 3 log a, (115) 
hy b uo 2 = Q a a~ 5 , h a b V b a 2 = E a a- 7 |loga| 2 .(116) 

Since A a — 0, £3-£~7 are undefined. So, to preserve lin- 
earity, all £i-£7 must be much less than 1. 



£1 



£3 



a X b\ _: 
Td ■ 



\nMZ a \ 
0, £4 



£2 



S a X b\ _ 



|2fi„| 1 

2 A a| 



|2E„| 



Ioga| (117) 
a"5|l oga |2 ( 118 ) 



£(, = 1 ' a 



1 FT" *' £l = lly I 

2 a I I 2 a I 



2 1 logal 



(119) 



So, in this case also some of the parameters £1, £2, £5, 
£6, £7 may become order 1, near the bounce. 



VII. CONCLUSION 



be adequate to give proper evolution of perturbations 
through the bounce. In sharp contrast with the result 
obtained in [l2j our result is independent of choice of 
gauge. We used gauge invariant variables that were not 
assumed to be small with respect to background. So one 
can evolve them through the bounce and match with cor- 
responding quantities in the expanding phase — but this 
would require the full nonlinear analysis. 

In this work, we consider only the contracting branch 
and used general relativity with usual matter distribu- 
tion as a correct theory to describe the dynamics of the 
universe. To investigate the nonlinearity of perturbations 
in a concrete manner, we have to take specific models of 
bounce. Currently we are working on this issue and hope 
to report in future. 
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Appendix A: Spatial Harmonics 



The tensor eigenfunctions (harmonics) of the spatial 
Laplacian V 2 = h ab V a h b c \7 c listed below fl5j].[29|. are so- 
lutions of the tensor Helmholtz equation. 



Previously it has been shown [ll|, [12( that in some 
commonly defined gauges linear perturbations grow so 
much near the bounce that that may invalidate the lin- 
ear perturbation theory. But the perturbations remain 
small if one uses some other gauge fixing condition. It has 
been argued that the commonly defined gauges are not 
well defined near and at the bounce. Hence linear per- 
turbation theory is valid if one uses only the well-defined 
gauge. Naturally, a question arises whether these results 
are real or gauge artifacts. 

In order to investigate this issue we used the covari- 
ant approach to perturbation theory. We focus on the 
evolution equations for density perturbation X a . The 
validity conditions of linear approximation of the (non- 
linear) density perturbation equations are set in terms of 
some linearity parameters. Then the linear perturbation 
equations are solved for a collapsing FLRW background 
near the bounce. The solutions are used to compute the 
linearity parameters. It is found that some of those pa- 
rameters grow beyond order unity near the bounce in 
both radiation and dust dominated cases. That means 
the nonlinear terms are comparable to the linear terms. 
So unless some special initial conditions are imposed on 
the variables such as shear and vorticity, perturbations 
may not be linear near the bounce. 

Thus we conclude that perturbations may not be linear 
near the bounce and linear perturbation theory may not 



V Qab...c 



0. 



(Al) 



(1) Scalar harmonics: Harmonics constructed from 
solutions of the scalar Helmholtz equation, 



V 2 Q(°> + ^Q(°) = 0. 



(A2) 



Vector and tensor eigenfunctions constructed from the 
scalars are 



Q(°) = -^> 6 Q(°), 



n (o) _ / a \ v n (o) 

= m 2 v <a v fc) Q(°> 



(A3) 



= (I) 2 h { c a h«V c V d Q^ + l -h ab Q^. (A4) 

(2) Vector Harmonics: Harmonics constructed from 
solutions of the vector Helmholtz equation, 

V 2 QW + ^QW=0, V»QW=0. (A5) 
Tensor eigenfunctions constructed from the vectors are 
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(3) Tensor Harmonics: Harmonics constructed from 
solutions of the tensor Helmholtz equation, 

Vq2> + £q<$ = 0, V b Q%=0, Q^ = 0. (A7) 

Q^-©V,4'. (A6) 
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